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SUMMARY 

The method of Poggi is used to calculate, for perfect 
fluids, the effect of compressibility upon the flow on the sur- 
face of an elliptic cylinder at zero angle of attack and with 
no circulation. The result is expressed in a closed form 
and represents a rigorous determination of the velocity of 
the fluid at the suiface of the obstacle insofar as the second 
approximation is concerned. 

Comparison is made with Hooker's treatment of the same 
problem according to the method of Janzen and Rayleigh 
and it is found that, for thick elliptic cylinders, the two 
methods agree very well. The labor of computation is, 
moreover, considerably reduced by the present solution. 

The third approximation to the compressible flow about 
circular cylinders, including the terms involving the factor 
(''o/coY, is also obtained and compared with the result given 
by Poggi. It is found that the expression given by Poggi 
is incomplete with regard to the terms containing the factor 

OoA-o) 4 . 

INTRODUCTION 

The purpose of this paper is to employ the method of 
Poggi (reference 1) to determine the effect of compressi- 
bility on the flow about elliptic cylinders. This prob- 
lem has already been considered by Hooker (reference 
2) who made use of the method of Janzen and Rayleigh 
but, owing to the necessity for expanding a certain 
function in the analysis, the "thickness ratio" of the 
ellipse to which his result applies is limited. The thick- 
ress ratio of an ellipse is defined as the ratio b/a, where 
u and b are the semimajor and semiminor axes, respec- 
tively. The method of Poggi, on the other hand, not 
only permits an unrestricted thickness ratio but also 
reduces the labor of computation. 

Briefly, it may be said that Poggi considers compres- 
sible flow to be replaced by an incompressible flow due 
to a distribution of sinks and sources throughout the 
region of flow. The strength of the distribution in the 
plane of the profde is given by 

and in the plane of the circle, into which the profile is 
mapped by a suitable con formal transformation, by 



where 

r, 6 are the polar coordinates of a point in the 
plane z(=x-\-iy) of the circle. 

R, 8 the radius of the circle into which the pro- 
file is mapped and the angular coordinate 
on this circle, respectively. 

X=-; tv=-^> V'^-Tze' * is the velocity 

potential of the flow. 
v, the magnitude of the velocity of the fluid 

in the plane of the profile, 
c, the magnitude of the local velocity of 

sound. 

Poggi then finds that the total velocity induced, at 
any point P(R, 5) of the circular boundary by the fore- 
going system of sinks and sources, is: 

bv 2 vg djr 

Ar=:r f „. Jf X * d /. , v ,, sin{d-5)d\dd (1) 

^tJoJo 2X cos (0 — 5) + V) 

Poggi's method of approximating the compressible flow 
of a perfect fluid is based on the assumption that the 
incompressible flow is a suitable first approximation 
and that therefore the values pertaining to that flow 
may be substituted for v T , v$, and zr in equation (I). 
The value of Av thus obtained then represents the 
effect due to compressibility and is to be added to the 
already known value for the velocity of the incom- 
pressible flow. That is, 

t'comp = Vincomp T At) (2) 

It is to be noted that, in equation (1), the local 
velocity of sound c is not a constant but is related to 
the velocity v of the fluid in the plane of the profile by 
means of Bernoulli's equation and the equation of 
state of the fluid. Thus, if the adiabatic equation of 
state is adopted, 

H^St'-s')] 

where c 0 , v 0 are the corresponding magnitudes in the 
undisturbed stream and 7=1.408 for air. 

In order to facilitate the solution of equation (1;. 
it has been the custom to replace c by <•„. This simpli- 

•Jt5 
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ficafion of the problem may he justified by the following 
argument. It has been tacitly understood that no- 
where in the fluid must the velocity- of the fluid exceed 
that of the local velocity of sound since the incom- 
pressible flow has already been assumed to be a good 
first approximation and the effect of compressibility 
is merely to distort the streamlines associated with the 
incompressible flow. As the maximum fluid velocit}- 
oecurs at the surface of the obstacle, there exists a 
value of r u -/c u 2 for which the maximum fluid velocity 
equals that of the local velocity of sound. -This 
critical velocity of the fluid is obtained from equation 
(3) by replacing v by c. Thus 



Co/' 



(4) 



This value for c is a lower limit under the condition that 
nowhere in the fluid is the local velocity of sound ex- 
ceeded. The maximum value of c occurs at the stag- 
nation point v— 0 and is given by 



(5) 



Tim- boili (lie niaxiinuni and the least values off occur 

on tin obstacle and everywhere else c nttX ^>c^>c lrttst . It 
follow fioin equation.- (4! and (.5) that 



r •> C I- 1 \ r . 

(',, 



= D.l)^s7 1+0.204 



Co' 



<0 



u.OdTo. The foregoing discussion thus shows that c/c 0 
may, »« a first approximation, be taken to be unity. 
Equation (1) then bocomer- 



l hi: 




r, 



Or _r» dr 



FLOW OF A PERFECT COMPRESSIBLE FLUID 
PAST AN ELLIPTIC" CYLINDER 



Let the jT j)lane be the plane of the ellipse and the z 
plane be the plane of the corresponding circle. Then it 
i- \sell known that Ihe Jmikowski transfonnat ion 



uliii'h increases very slowly as r 0 ,'c 0 approaches unity. 
In fa« i. it is seen that the upper limit for c "" r ~ f ^ 



which it. is transformed are, respectively, and 



i?- 



a- 

Tt' 



The thickness ratio t then becomes: 
B- " 



f= 



Ii l-< 



R - r R 



1 -f CT 



or 



1-/ 

T-H 



where 



If w denotes the complex potential of the incom- 
pressible flow in the z plane when a stream of velocity ^ 0 
impinges on a circle of radius R in the direction of the 
negative x axis, then 



(8) 



The complex velocity in the f plane is then given by 

dw_dv:dz 

rff-rfz d{ 

or 

rfir_ c 2 -i?- 



(9) 



When X — — and cr = ( ,' ) are introduced, it follows that 



tor? 



To' 



1— 2X : cos 20-j-X 4 
l-2<rX 2 cos 2e-r-tr 4 X 4 



(10) 



Following Poggis procedure, the Fourier development 
of tr/vfvn)} be obtained. Thus, by the use of the 
expansion 

1 

1— 2<r'X 2 cos 20-fVA 

(see appendix, sec. 1), 
it follows that 



i- i 

—i=o a o+^Zo 2n cos 2nd 

<0 - ri=i 



(in 



where 



1-f (l-2^ )X j 
1 — o<\ 4 



7^ I 



(7) 



map- • ! t uric ol indiii- " with its center at the origin 
of the : j)lane into the line segment ( — 2a. 0; 2n, <),) in i 
I In- plain . .\l-o. the circles concentric w ith the circle i 
of radio- a are transformed into a family of confocal i 
ellip-ev with common foci at (.— 2a, 0) and {2a, 0). Jf j 
/i't <;,i denotes the radius of one of these circles, then ' 
the semimajor and semiminor axes of the ellipse into | 



and for i<= 1 , 2, 



«-- = 3 V ,. V1 M _^ <t (a\y 



a 2 (l — <7 4 Vj 



Also from equation (S) 



r,= — r 0 (l — X 2 ) cos 8\ 
re=t; 0 (l -f- X 2 ) sin 0 | 



(12) 
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Then, inserting the expression? for >\, and Ci given 
by equations ill) and j 12 1 into equation H)'i and making 
use of the integrals 



Jo l--2\t 



sin <n — ")i 

. —.Mil r'hlA = 

OS . '/ - 0 I — A" 



r < 

[_jrX 



0 if /i=0 

"~ l ens no if n > I 



-cc appendix, -e<\ I [ i 

j'-' T <in ■>)-,]■ r (i if n =0 

.1 1^2X ..;-5. -,V '■' w "'^-l_-.T\"--in /,3if«£l 

it follows without diliieulty that 

A ' ; =^r -in o|S f 1 ) sin <2« + l>5 
in - L '< -" 

| 1 (,V" r, f/..,-X-"-'« Jn .,i</\j (13) 

w lie re 



(1 + <r) (1 — <rr . 1 -J- 2 or cos o — <r 

- £ sin 25 log 5 > 

2cr 1 — Jo- cos d — a 

, 1~~ Q~Y . i i , 2<r sin o 

+ — — U — c 4 ) cos 25— 2cr tan 1 — _~ -.- 

+ 2^(1 — o-'-'-rT 5 ! ^in 5 — c 2 sin 35 Jj ) 



1 *J 



For 5= .7. the position of maximum velocity on the -ur- 

l 

; face of the elliptic cylinder. 

t 

/ At\ ml — <rV I— <r pi — <ri-", I - ,r J 



Substituting for the a ln 's from equation (11), equation 

(V.]) takes the form 



1)5 



~~$ l —^f^n 5_(i_^2(o fl j.i) sin (2/i+ 

j:: h 

Replacing X-' by r. for purpn-cs of integration only, it 

follows that 



/=S (2»~- 1 ) *in (2/i+ ho ( '! ~^*-! U 2 X') VX 2 

n=0 ' Jo l— o- A 

= -/?./>. of /SC-'/i-rl V-''- 1 ' 5 ( '{"-pfan^r 
«=o ' Jo 1 — o-V 

Jo (1 — o- r-«- ,J )- 1— <rV 

= ,, - 1 ... ; 4 ,.. | (1 rT" )2 [(l-:3+-o- < ) sin 5 
(1 — 2o~ cos 25 + o- 4 }-( 2c- 

+ <r sin 35] log ~p; 

i 1 — °"r/, i i\ o «i l ,2a sin 8 
-\ [(l + o- ) cos 25 — 2c"] tan 1 — >- 

0" 1 — o" 



! — rr- 



tan 



or 



(1 +«*)(! -a 3 ) 



l+2<r cos 54-^ 



s sin 2o log : — ^ —, — , 

2tr 3 1 —2a- cos o + cr- 

+2[(l + cr+o- 4 ) sin 5-c 2 sin 35]j 

Therefore 

&v ill — a 2 / . . 
- = .-, — -, — sin o 
in 2 <r" V. 

— 7,— rr - > ■-— .vTi — r» o ■> ~ (l+3o-+or 4 ) sin o 
(1 — 2o~ cos 2o + <r'j-i 2ff- (_ 

+ o- 2 sin 35 losj ! — — , 
i ■: , .i77>— :;u —it 



It is interesting to note that the expression for 2\r r„ 
at the surface of a circular cylinder fixed in a stream of 
velocity c a impinging on it in the direction of the 
negative x. axis may be obtained from equation i'l4) by 

allowing <^(j=^) to approach zero. Thus, making use 
of the expansions 

( 1 -2o- 2 cos 25+ o-V = < \^4\ 1 * ^ + M ' 
— (/i— l)o- 4 Jo- 2 " cos 2/'5- t<ee appendix, see. Ill) 

. l+2o- cos 5 + c 2 a'"' 1 , , . 

log — 5 — — i-42J 0 . i c05 l-"-rl "> 

° 1— 2<j cos o + c- ,7^/2/! + 1 

2a sin 5 



tan 



.. sa sin o .^-v <r~ ' . , , . 

1 -j - = 2>_,v — — ran (,2/i+Uo 

1 — 0~ n =0-« T" I 



log r, = 2> ,^ - 

3 l — o- ^o2/i.— 



it follows, neglecting terms containing powers of c 
higher than the second, that 

~ = ^{l -^ 2 )(| sin 5-sin 3o) 



or 



Li m-^=y|sin5-isin3o") fit)) 

This expression for Av/v Q agrees with that obtained 
by "the methods of Janzen, Rayleigh, and Poggi (ref- 
erence 3). 

The effect of compressibility, i. e., Ai'/tV), having been 
found, it follows according to equation (2) that the 
total velocity at the circular boundary in the z plane 
is given bv 

(J-) =2sino + f (17) 
and on the elliptic profile in the j" plane by 

(,^L/.p« = (I -2a- cos 25 4- c'j '-C T o ) ttrlh (l8) 
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Table ] shows the comparison between the values of 
^V' , o)<ii(7-jc calculated according to equation (18; and 
t J lose obtained by Hooker for an ellipse of thickness 
ratio t = ]'i or o'-~){. The values for the corresponding 
incompressible flow are included. It is seen that the 
results of (he two methods agree very well. This agree- 
ment is not unexpected since Hooker's method is par- 
ticularly applicable to thick ellipses. Consider, how- 
ever, a slender ellipse, say <=Ko or <r=Y u . Table II 
shows the comparison between the exact calculations of 
the present method and the results obtained according 
to Hooker's method. The disagreement is more evident 
than that shown in table I for the thicker ellipse. 
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THE PRESSURE DISTRIBUTION 

According to Bernoulli's theorem and the adiabatic 
equation of state, if p and p are the pressure and density 
of the fluid, then 

where p 0 and p 0 are the pressure and density, respec- 
tively, in the undisturbed stream. Expanding the 
right-hand side, of the foregoing equation and neglecting 




IOO 



Figure 2.— The pressure of the fluid on the surface of an elliptic cylinder of thickness 
ratio 1/10 for compressible and incompressible flows withjo/co=0>" 
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forms invol viim powers of n higher than the first yields 

FM-$M-&)'+ • ■ • • <•»> 



The procure distribution is rhen obtained by sub- 
-tituting for '■ •■. fru'-i ■■qnar'nm . !s>. Table III shows 
the pressure di->:ril>uiiuii o\<-r tin* -urlace of an ellipse 
i>f thickness ratio I 10 with ' <■. -0 <~u , and figure 2 
allows the graph of this distribution together with the 
one due to the corresponding incompressible flow. 
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THE ATTAINMENT OF THE LOCAL VELOCITY OF SOUND 
AT THE SURFACE OF AN ELLIPTIC CYLINDER 

According to equation (4) the critical velocity of the 
fluid is given by 

[r^)- y Ti:u^y+i (20) 

For an elliptic cylinder, at zero angle of attack, the 

critical velocity occurs at 5=^; the position of maximum 

velocity on the cylinder and also in the region of flow. 
Hence substituting from equation (IS) for (e/» 0 )«iifp« at 

o=^ yields a cubic equation in the variable m- 
Thus, from equation (15), if 



then 



0 



_9- 



-tan-'ff+S 



f/(cr)]V+4/(^M 2 d-[4-^ r j(H- ( r 2 ) 2 ], i 
2(1- 



= 0 



(21) 



where 7=1-408 for air. 

Table IV gives the critical values of v 0 fc 0 for the 
entire range of thickness ratios including the limiting 
cases of the straight-line segment and the circular 
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profile. Figure 3 shows the critical values of >'rjc n ( — \ a) 
plotted against the thickness ratio. 



^-Straight-line segment 
I I 




Thickness r 



Figi-re 3.— The critical r.\tio r./fo as ;i fusu-ti .n .if tbe tli:ck:i«;s rjt o (. 
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THE THIRD APPROXIMATION TO THE COMPRESSIBLE 
FLOW ABOUT CIRCULAR CYLINDERS 

In reference 2, the opinion is expressed by Hooker 
that the terms involving (r 0 /c 0 )\ thus far neglected, 
may become of considerable importance as the local 
velocity of sound is approached on the ellipse. Hooker, 
however, did not investigate the matter any further. 
In reference 4, Poggi calculated these terms for the 
compressible flow about a circular cylinder, but a close 
examination of his work shows that not all such terms 
were taken into account. In what follows the terms 
neglected by Poggi will be obtained and compared with 
the alreadv existing ones. 
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The fundamental integral equation (1 I may ho written 
folk Ws: 




r 0 d\ 



1 'lj Tu~ 

X F 0 ~du 



1-2X cos (0— 5; + X s 

[l -^^(l ~^»)+ • • ] si » (e-t)dM6 (22) 

where l,'c : ' has l>een replaced by a power series in 
/i ( — i\r. : c { f) obtained from equation (3); i. e.: 

H,D-^0-S>---] 

The method followed by Poggi was to substitute for 
r M r„, and r expressions pertaining to the incompressible 
How and thus obtain the following result: 1 

7"!, - (J sin 5- r, sin 35^+ (o - 1)( ^ sin 5- ~ sin 35 



(23) 



Tlic velocity for the compressible How at t he surface of 

the menhir cylinder then becomes: 



J \ntou,p ) 



Ac 



(24) 



win n A> /•„ i- ghen by equation (23). 

h<|ii:ition >2i) i hut- represents the second approxima- 
iioi' in the <•( >mpre->ible flow, the first approximation 
Ix n><.' tin puich incompressible How given by r lneom ,,jv n . 

The third :i])proximation may be obtained, at least in 
principle, by >-iilvMituiing for r„ <>, and r 2 in equation 
12:' i cxpre— »i««n- based <m the second approximation. 
Such expre-- ions, ii< far as the terms involving are 
concerned, are given in reference 3 and are as follows: 



r, 



(1-X-) cos 0-mU — |4 x *+I v — ^ x ") cos 6 
+ (;>-jx«)«*:w]-^.. 
fl ■ X-» <in 0 - pjj |o x2- ^ x4 -^T^ x6 ) snl e 

-(- : ;x--;r)>i„3^.. 

- ( - % >• ,\ .!x v )<-o-2« ; X 2 cos 4«Tj- . . 

, • ; mi - .fl -l ^ w ere obtained for elhpnc cylimler- 

r<"iii' • r- '*'•-< :m equation <2i'i for the circle llou- 

in. tin urn- n! ',:iij:i<: -in34are sivoa, respectively, as 

. -l:s.iii i rror i:> ( iK ul iSiuu- ]'r>lole.<: (refercnce3) pives for the 
i ii 'I ci. it of mi i ii.i i.ilu' ' .i «ln • . in we.', of the indepeaderit check of Posprs 
r» *j 1 1 t>\ ih( ifin, r . i- i>che\«'.i tn U .. ini-i«r::*T 



a'i'i 



f. ,!• 1 



where it is recalled that X - 

/■ 

When the foregoing expressions are substituted into 
equation (22) and only the terms involving n and /r are 
evaluated, it is found that, besides the terms given In- 
equation (23), the following ones involving must be 
included: 

These terms seem to have been overlooked by both 
Poggi and Pisfolesi (reference 3). 

The third approximation to the compressible flow at 
the surface of the circular cylinder then becomes: 



To 



5 + (j sin 5--|sin 35].u 

, T/37 . 25 . _ , 3 . . \ 
\40 sin "M sl " "*"8 SU1 / 

/ 23 . 11. 



+ g sin f>5 



(2u) 



It is interesting to compare the magnitudes of the 
various terms in equation (2ii) at the position of maxi- 
mum velocity 5=tt/2 and for the critical value fi=0.1670 
(obtained b\ means of equations (20) and (26)). Thus 



(2 sin 5) , = 2 



sin S-^ sin 35^_. = 0.194S 
sin sin 36 + ^ sin 55 ) a _ x = 0.0fi53 

■(7-l)(p^ 0 sin S~ sin 35 + -J sin :,5 ). r =0.00C7 



Thus, it is seen that the terms in\ olving ^- do become 
of importance with regard to the p tenn> as the local 
velocity of sound is approached on the circle and that 
the main contribution i- made by f\pr<— sion t2")). 
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APPENDIX 

! Therefore 



1 — -J-r-.V cn~ 2'J— o- ! X ! 
I f 2 co- "J" I- - i!:ti . I liv t - " '. i il-'li 

I ! 



11 The Integra W 



1— 2<rX J ■ -1 — r.\ . ' ; l-ff-.W'-") 

Since, by the biiinmi.il theoiein, 

and 

(I— ff-'W-"!- 1 -^f.r-'.V. » 

..- -ii 

it follows tint t 



-in >6—o) 



and 



! *-J."i=35 



„ I— 2X cos i# — oi - X 
-in (d—o) 



COS /(#'/<' 



, sin n6<lt) 



Lot 

a nil therefore 



77= S S^^'Vy-"' 

j=0 fc=0 



j+Jfc = n 



j—L^n—2h; j=n—k 
The double series tlien boeomc*> 

//-S v\, T - X vv- ■'- - 
« -ii ■. =ii 

Tlie terms of iliis -erics can he grouped in pairs such 
that 

n .i - 1 



!X cos (0— S)-~ \ 2 
If 2cos (9—5) is replaced by e ,w_< ' fe -, ^" M , 



then 



1 



1 



I — 2X cos (d—o) -f- X- > l — \e 11 ' 



I — X«~ 



(1) 



e u«-a) _ e -..9-s)| 1 _ xg'«-«) i — Xe-HO-si) 

2i sin (0— o)|_;^o ' T^/ 
Therefore 

./,-;-,•,/,= 1 . 1 "fe x-v," "-«• x • 

-?J0 -1-0 r.=0 

lie placing r* by :, 

- J =o 



, «/}—[.. .... ,. 

where ^ or — ~ } — i.- the upper lmut according as /( is 

even or odil and where the factor 2 is omitted from the 

term for which ft is even and This term is 

independent of 0 and there is only one such term, not 
two. 

The scries (I) may be written as 

n a = l 

77= S S An cos (»-2Jb)20 

li =0 t: =0 

where 

.4„=2(<r 2 X 2 )" 

Expanding this series and rearranging the terms in the 
form of a Fourier series, 



>ince 



j 2-i, when p— — 1 



But 



and 



tf=^S^U,-fZ; cos anffS-^M. 



1 = * 1 

iS--U = S(<r 4 X^=-^ r 
2„=o n=n 1— c X* 



. : .--2(VX-|- , ^Ct'X 1 )* = 



it follows that w = /< — 1 and therefore 

<-,/,=-/ x n -V'" 5 

Hence, for 1, 

</,= — -X' !_! sin «<5 and ./, = -X" _1 cos //<5 

III. The Fourier Expansion of 

1 



(l-2o-X i cos 20-r-o- 4 X 4 ) 2 
In analogy to section I, replace 2 cos 22 by V 
Then 

1 



7/ 2 = 



(!-2<rX 2 cos 2d+o- i \ i y 
1 



( I - o-'-'X-e-'"; -'( L — <r W-'" ) 2 
According to the binomial theorem 



and 



1 1 - <7 2 XV">) == 2(jJ- 1) (<r* \W" 



( 1 - .r-'X-e --■") ~ l = SfH 1) (oW) V J,t * 
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Therefore But 



Let 

and therefore 
Then 



j+k=n 
j—k=n — 2k, j=n—k 



H 2 =J2 J}(n-k+l ) (Je+1) 0 2 X 2 ) V"-** 1 ' 

71=0 k=0 

The exponent of c is 2i[(n— k)— k]6. If k and n— fc are 
interchanged, the exponent of e changes sign but the 
coefficient of e remains unaltered. The terms can 
therefore be grouped in pairs so that: 

n v_— ] 

ir-^-lit (Ti-Ar-i-lXH-lK^X 2 )" cos (n-2i)2fl (4) 

where the factor 2 i~ omitted from the term for which 

i, \< c\ en and /. — ~ 

Tlie -eric- (-J; may be written a;- 
t. )■ — i 



V hel i' 



]]'^J2 A„.,. cos (n-2k)26 

n 0 /. = u 



A. ,=l'0'-/.-+l)(A-+l)(ff 2 X s )" 

K\|)iiii(IiiiL r thi< series and rearranging the terms in the 
t'li'in of ;i l-'uurier serie-. 

/A - \1b >.-r S cos 27,0 S * 

-J-0 71=1 A=0 



and 



1 » » 

oS-^2t.t=S(^-rl) 2 (^ 4 X'')*= 7- — 1 — 4 j 

^t = 0 *~0 (1 — (TAJ 



S^^.t=2(,7 2 X 2 )' 1 2:(n+Z:+l)(t-Ll)( C r 1 X 4 )' 

t-0 *-0 



— 2(<r 2 X 2 ) - f ,- . T , 



Therefore 



# 2 = 



1 



(i-F^x 4 )-f2S[(«+i) 

71 = 1 



(l-<r 4 X 4 ) 

-fa-lVX'JtVX 2 )" cos 2?,e 
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